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SPATIAL ASSOCIATION BETWEEN A SPIDER WASP AND 
ITS HOST IN FRAGMENTED DUNE HABITATS 
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ABSTRACT, In patchily distributed habitats, species potentially occur wherever conditions are suitable 
or show a restricted distribution, influenced by patch quality, geometry and configuration. If patch isolation 
appears to be the main determinant of the species’ distribution then dispersal ability is supposed to be 
limited. Although only scarce literature is available, dispersal limitation seems to be an important factor 
in determining the spatial population structure in spiders. In this paper, we document on the spatial 
population structure of the rare wolf spider Alopecosa fahrilis , restricted to fragmented grey dunes along 
the Flemish coast (Belgium) and ask whether its distribution appears to be affected by aspects of patch 
configuration. Simultaneously, we investigated whether the local distribution of its main parasitoid, the 
spider wasp Arachnospila rufa (Hymenoptera, Pompilidae) was associated with its host. Our results in¬ 
dicate that A. fabrilis shows an aggregative population structure, which is determined by the distance to 
nearest occupied patch, indicating that spatially correlated habitat quality probably determine its occur¬ 
rence. Although spider wasps are generally characterized as non-specialists, the almost complete covari¬ 
ation between its spatial occurrence and that of A. fabrilis , indicates that spider hunting wasps may, at 
least temporally and locally, show a restricted host-range. As a result, the presence of a rather generalist 
parasitoid is a good predictor for the presence of nocturnal and burrowing dune wolf spider. 
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In fragmented habitats, species live either 
in patchy populations or in metapopulations 
(Harrison 1991; Hanski 1999). In patchy pop¬ 
ulations, individuals move freely among hab¬ 
itat patches, while in metapopulations, most 
individuals stay in a single patch during their 
entire life, but dispersing individuals enhance 
strong colonization-extinction dynamics. This 
results in a population structure in which some 
suitable patches remain vacant. Besides the 
colonization abilities, changes in habitat qual¬ 
ity can also attribute to local extinction dy¬ 
namics, as demonstrated for specialized but¬ 
terflies (e.g. Thomas et ah 1992; Ravenscroft 
1994; Moilanen & Hanski 1998; Bergman 
1999) and backswimmers of the genus Noto- 
necta (Briers & Warren 2000). Therefore, 
studying the spatial structure of populations 
during successive years by considering as¬ 
pects of patch geometry, quality and config¬ 
uration enables us to assess actual dispersal 
limitation and population dynamics in an in¬ 
direct way. Snapshots of patch-occupancy in¬ 
cidences provide an alternative indirect meth¬ 


odology to study dispersal abilities or 
populations viability with respect to popula¬ 
tion size (patch size) resulting from historical 
habitat fragmentation and population dynam¬ 
ics. As shown for a coastal dune wolf spider, 
living in fragmented grassland habitats, patch 
occupancy patterns may depend both on as¬ 
pects of habitat quality and on different modes 
of dispersal through the surrounding matrix 
(Bonte et ah 2003a). For spiders in general, 
cursorial dispersal and ballooning induce dif¬ 
ferent colonization and extinction patterns, de¬ 
pendent on the species’ mobility, its niche 
breadth and Its propensity for aerial dispersal 
(Bonte et ah 2003a, 2004b). Especially bal¬ 
looning dispersal appears to be important for 
occupancy patterns at extended temporal and 
spatial scales (Bonte et al. 2003a, 2004a), 
while cursorial dispersal induces short-term 
colonization events at small spatial and tem¬ 
poral scales (Bonte et ah 2003a). For spiders, 
direct estimates of realized dispersal are rare 
and indicate dispersal distances up to few- 
hundred meters within suitable habitat (Krei- 
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ter & Wise 2001; Bonte et aL 2003a; Samu et 
al. 2003) but restricted dispersal up to few me¬ 
ters in unsuitable, often densely vegetated 
habitat (Bonte et al. 2003a). Indirect estimates 
of dispersal, estimated by population genetic 
analyses, are more common (e.g. Boulton et 
ai. 1998; Ramirez & Fandino 1996; Ramirez 
& Haakonsen 1999; Gurdebeke et al. 2000). 
In contrast to the latter, studies on occupancy 
incidences, ideally conducted during several 
successive years, reveal indirect estimates of 
dispersal that are independent of historical 
bottlenecks and selection pressures in spatially 
structured habitats (Peterson et al. 2001; Bon¬ 
te et al. 2003a). 

In addition to earlier studies on the mobile 
wolf spider Pardosa monticola (Clerck 1757), 
we here report on the population structure in 
coastal dune habitats of the specialized fos- 
sorial wolf spider Alopecosa fahrilis (Clerck 
1757) with limited aerial dispersal abilities 
(Bonte et al. 2003b) and of its main (winged) 
parasitoid, Arachnospila rufa (Haupt 1927) 
(Hymenoptera, Pompilidae), knowing to host 
on larger Alopecosa wolf spiders (Koomen & 
Peeters 1993; Peeters et al. 1994). According 
to theoretical work, host-parasitoid metapo¬ 
pulations dynamics strongly interfere and oc¬ 
cupancy patterns should strongly overlap in 
case of limited host interpatch dispersal abil¬ 
ities and relatively low infection rates and par¬ 
asitoid survival rates (Hassel et al. 1991; 
Comins et al. 1992; White et al. 1996; Hanski 
1999). 

The population structure of parasitoid and 
host were investigated in the Flemish coastal 
dunes, where grey dune habitats, being the op¬ 
timal habitat for both A. fabrilis and A. rufa , 
are truly fragmented since the Second World 
War (see e.g. Bonte et al. 2003a; Provoost & 
Bonte 2004). This was done in order to test 
the following hypotheses: (i) A. fabrilis has 
limited dispersal abilities, resulting in an iso¬ 
lation-dependent population structure, (ii) A. 
fabrilis , having low population densities only 
occupies larger habitat patches and (iii) local 
and temporal host-parasitoid association result 
in similar distribution patterns, although Pom¬ 
pilidae have better dispersal abilities and are 
believed to be rather generalistic in prey 
choice according to prey species (Finch 1997). 
Additionally, as A. fabrilis is nocturnal while 
its parasitoid is diurnally active, we asked 
whether the presence of a day-active parasit¬ 


oid could be an indicator for the presence of 
populations of a nocturnal, fossorial spider. 

METHODS 

Study Area, Study Species. —Fieldwork 

was conducted in the Flemish coastal dunes, 
located between the cities of De Panne (Bel¬ 
gium) and Bray-Dunes, France (51°05'N, 
2°32'E) consisting of 52 discrete grey dune 
patches, varying between 0.05 and 27.6 ha 
(Fig. 1). Grey dunes, known as “Fixed coastal 
dunes” are most readily defined using plant 
communities. Vegetation includes Atlantic 
moss dominated dunes (mainly Tortula rur- 
alis) as well as dune grassland (with a distinct 
organic soil layer) belonging to the Cladonio- 
Koelerietalia syntaxon in case of lime rich 
grey dune, and to the Trifolio-Festucetalia ovi- 
nae syntaxon in case of decalcified grey dunes 
(Provoost et al. 2002). In this study, only grey 
dunes without substantial soil development 
were surveyed since they are the habitat of 
both study species (see further). Ecologically 
it is merely the dry component of the 
“stressed dune landscape”. The main differ¬ 
entiating processes are related to dune fixa¬ 
tion, soil formation and vegetation develop¬ 
ment (Provoost et al. 2002). At present, rough 
grass and scrub encroachment result in a se¬ 
vere fragmentation within a matrix of dense 
dune vegetation (shrubs, dense grassland). In 
an earlier paper, we identified typical spider 
species for this habitat (Bonte et al. 2002). 
One of the most habitat-specific species is Al¬ 
opecosa fabrilis, although its overall Indicator 
Value was low, hence, indicating a relative 
low occupancy rate. 

Alopecosa fabrilis is used as model organ¬ 
ism in this study. It is the largest lycosid spe¬ 
cies, living in self-made burrows in dry sandy 
habitats from temperate regions in Europe 
(males: 10-12 mm; females: 13-16 mm; Rob¬ 
erts 1998). The species has a nocturnal life 
style in which especially males leave their 
burrow in search for mating partners during 
autumn, with a peak activity during Septem¬ 
ber (Roberts 1998; Bonte, unpub. data.). Fe¬ 
males are active during autumn, but especially 
during late winter and early spring. The spe¬ 
cies occurs in low densities, has a biannual 
life cycle (Bonte et al. unpub. data) and does 
not perform ballooning dispersal under labo¬ 
ratory conditions (Bonte et al. 2003b). 

During its period of activity, the spider 
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Figure 1.—Map of the surveyed grey dune habitat patches in the coastal dunes of De Panne-Bray 
Dunes, Belgium. Black: patches occupied by Alopecosa fabrilis ; white: unoccupied patches. Circles in¬ 
dicate patches not occupied by Arachnospila rufa. 


hunting wasp (Hymenoptera, Pompilidae) Ar¬ 
achnospila rufa has been observed to be the 
main predator (Koomen & Peeters 1993; Pee- 
ters et al. 2004; Bonte, pers. obs.; Nieuwen- 
huijsen, pers. comm.). Alopecosa schmidti 
(Hahn 1834) has been recorded as prey in 
Middle-Europe (Schljachtenok 1996). Female 
pompilids provide each cell of their nest be¬ 
low the sand surface with only a single par¬ 
alyzed spider, on which one egg is laid. Ar¬ 
achnospila rufa is large (body size up to 18 
mm) and common in sandy regions of Bel¬ 
gium and the Netherlands. It reaches adult¬ 
hood from June-October (Peeters et al. 2004; 
Nieuwenhuijsen in press). More detailed in¬ 
formation about its ecology and life history is 
unfortunately not available. Voucher speci¬ 
mens of both Alopecosa fabrilis and Arach¬ 
nospila rufa are deposited at the Royal Bel¬ 
gian Institute for Natural Sciences. 

Field Survey & Statistical Analyses.—Oc¬ 
cupancy patterns of A. fabrilis were recorded 
from 25 August-7 October 2003 using pitfall 
traps (diameter 9 cm, 6% formaldehyde-de¬ 


tergent solution). In each of the 52 grey dune 
patches, at least five traps were randomly 
placed, depending on the patch area. Patches 
were digitized from aerial orthophotographs 
with a Geographic Information System 
(Arcview 3.1) and discrimination of vegeta¬ 
tion types was based on vegetation-specific 
red (RED) and near-infrared (NIR) reflectance 
values (Provoost et al. 2002). Because of its 
high activity during this season, the use of pit- 
fall traps have been shown to be very useful 
in catching this typical spider species living 
within this habitat (Bonte et al. 2004c). The 
presence of the hunting wasp was recorded in 
the same pitfall traps, but completed with de¬ 
tailed field surveys during sunny days. Al¬ 
though not all observed specimens were col¬ 
lected to confirm identification, individuals 
could be identified in the field by their large 
size and distinct abdominal coloration. From 
GIS, we measured patch area, the distances to 
the nearest suitable patch, as a measure of 
patch isolation and distance to the nearest oc¬ 
cupied patch (nearest-neighbor distance) as 
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measurement of population isolation. For A. 
rufa , we did not use the latter isolation mea¬ 
surement because its presence was only re¬ 
corded during September, hence ignoring pos¬ 
sible different distribution patterns during the 
previous summer-period. 

Patch occupancy patterns for A. fabrilis (0: 
vacant; 1: occupied) were analyzed by logistic 
models for binomial data and logit-link with 
backward elimination of the non-significant 
parameters (SAS 9.2). Patch area, patch iso¬ 
lation and the interaction between both were 
included as independent variables. For A. rufa , 
the occupancy status of the patch by A. fa¬ 
brilis was used as an additional categorical 
variable. 

RESULTS 

Patch occupancy of Alopecosa fabrilis was 
not significantly influenced by patch area (x? 
= 2.255; P = 0.133). The distance to the near¬ 
est suitable habitat (x 2 = 1.168; P = 0.280) 
and the interaction between both (xf = 0.208; 
P = 0.648) did not contribute significantly. If 
the nearest-neighbor distance was taken into 
account, a significant negative relationship 
was found between isolation and patch occu¬ 
pation (estimated slope:—0.035 ± 0.014; x? = 
6.322; P = 0.012; Fig. 2). Other parameters 
remained non-significant. Although the model 
describes patch occupancy patterns in a sig¬ 
nificant way, 32.7% of the cases were mis- 
classified. Mainly patches near the inner dune 
front were occupied (Fig. 1). If only patches 
at the inner dune front were used for analysis, 
neither distance to the nearest patch (xf = 
0.028; P = 0.865), distance to the nearest oc¬ 
cupied patch (xf = 1.144; P = 0.285) or area 
(Xi = 0.476; P = 0.490) explained distribution 
patterns. 

After backward elimination of non-signifi¬ 
cant parameters, the occupancy of Arachnos- 
pila rufa was only determined by the presence 
of Alopecosa fabrilis (xi = 34.41; P < 
0.0001). No other lycosid species were re¬ 
corded during the survey period. Landscape 
geometry appears to be unimportant during 
this season since the explanatory power of 
both patch area and patch isolation are ex¬ 
tremely low (all variables and interactions: x 2 
< 1.29; P > 0.256). Patches occupied by A. 
fabrilis but not by A. rufa (n = 4) during this 
period were significantly smaller than those (n 
= 11) occupied by both species (Mann-Whit¬ 


ney U-test: Z = 2.74; P = 0.006; Fig. 2). Ar¬ 
achnos pila rufa was not recorded in patches 
without A. fabrilis populations. 

DISCUSSION 

Population structure of Alopecosa fabril¬ 
is .—Our results demonstrate that the distance 
to the nearest occupied patch contributed sig¬ 
nificantly to the patch-occupancy model, 
whereas the overall habitat structure (i.e. ge¬ 
ometry and configuration of all available 
patches) did not. As a result, the population 
structure seems not to be affected by local po¬ 
tential population sizes (related to local pop¬ 
ulation viability; Hanski 1999) and distant dis¬ 
persal abilities by ballooning. The latter is in 
agreement with the absence of aerial dispersal 
under laboratory conditions (Bonte et al. 
2003b). As only the nearest-neighbor distance 
contributed significantly to the patch-occupan¬ 
cy model, the population structure of A. fa¬ 
brilis seems to be clearly aggregative, indi¬ 
cating low abilities to disperse (and penetrate) 
cursorially through the hostile matrix or a pre¬ 
dominant role of habitat quality. Our data, 
however, only result from a short term survey. 
But preliminary data on patch occupancy pat¬ 
terns from 2004 show similar distribution pat¬ 
terns (Bonte et al. unpub. results), indicating 
the absence of strong turn-over dynamics at 
extended, though still rather short time spans. 
The fact that patch area (related to local pop¬ 
ulation size) does not explain occupancy 
patterns seems to indicate that the spatial pop¬ 
ulation structure is more likely to be deter¬ 
mined by random extinction-colonization dy¬ 
namics within aggregative clusters of closely 
connected suitable habitat remnants, which 
appear to be situated near the inner dune front. 
As the interaction between patch size and 
patch isolation does not explain occupancy 
patterns, a rescue effect resulting from source- 
sink dynamics is unlikely (Dias & Blondel 
1996). 

The aggregated pattern may result from dis¬ 
persal limitation after historical population dy¬ 
namics or from spatially correlated habitat 
characteristics, related to habitat quality, as 
experienced by the species. Although both are 
possible, the second seems to be more prob¬ 
able as especially patches near the inner dune 
front were occupied. In coastal dunes, char¬ 
acterized by a successive ontogenesis, grey 
dunes, although structurally similar may show 
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Figure 2.—Relation between patch area, nearest-neighbor distance and occupancy incidence for A. 
fabrilis populations. Arrows indicate patches occupied by A . fabrilis , but not by its predator A. rufa. 


extensive spatial correlation in environmental 
characteristics related to microclimate, soil 
properties and aeolic dynamics (sand over- 
blowing). The fact that considerable non-spa- 
tial correlated variation in habitat quality af¬ 
fects spatial structural patterns has already 
been demonstrated by e.g. Briers & Warren 
(2000) and Bonte et al. (2003a). For a fosso- 
rial species like A. fabrilis, it is not impossible 
that e.g., grain size properties, soil water re¬ 
tention, soil formation and the intensity of 
sand overblowing strongly determine burrow¬ 
ing ability and the sustainability of created 
burrows. In the same coastal dune landscape, 
the grassland inhabiting wolf spider P. mon~ 
ticola shows a habitat-quality dependent pop¬ 
ulation structure in which both patch area and 
patch connectivity determined occupancy pat¬ 
terns (Bonte et al. 2003a), although not in an 
aggregative way. As vegetation structure is 
alike In all investigated grey dune patches, 
similarity in habitat quality for A. fabrilis is 
probably more related to pedological (possibly 
soil stability or differences in grain size near 
the inner dune front) or microclimatological 
(less buffered temperatures far from the sea) 


conditions and, possibly, the reason for the 
observed clustered occurrence. Also, the ab¬ 
sence of A. fabrilis in large patches more close 
to the sea, in which large populations should 
be able to persist, provides more likely evi¬ 
dence for a reaction towards spatially corre¬ 
lated environmental properties and optimal 
habitat quality in older grey dunes, situated 
near the inner dune front. In this case, dis¬ 
persal limitation is of inferior importance, as 
illustrated by the model in which only patches 
near the inner dune front were included. 

Host-par asitoid association.—Spatially 
explicit host-parasitoid metapopulation mod¬ 
els generate spatially correlated abundance 
and, hence, occupancy patterns In host and 
parasitoid (Hanski, 1999). Arachnopsila rufa 
and related spider hunting wasps, however, 
appear to be more generalist parasitoids on 
(larger) lycosid species (Koomen & Peeters 
1993; Endo & Endo 1994; Finch 1997; Pee¬ 
ters et al. 2004). The almost complete match 
between occupancy incidences of host and 
parasitoid indicate that, at least in the Flemish 
coastal dunes and during the survey period, A. 
rufa behaves as a specialized parasitoid on A. 
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fahrilis . This indicates that at least one species 
of pompilid wasp may behave as a host-spe¬ 
cialist during some periods of its adult life. 
According to White et ah (1996), spatial pop¬ 
ulation patterns in which both host and para- 
sitoid occur in the same patches are only static 
in case of low host dispersal and low parasit- 
old survival. Certainly the first assumption 
seems to be fulfilled in our study. However, 
A . rufa is already adult earlier in the season, 
and should, as a result, behave more flexibly 
in prey selection. As predicted by optimal for¬ 
aging theory, prey size selection should be op¬ 
timized in relation to the species 5 load-carry¬ 
ing capacity (Evans 1962; Coelho & Ladage 
1999), The latter authors indeed found that 
large female wasps searched for prey that 
matched their own body mass and lift capac¬ 
ity, As a result, A. rufa should prey on simi¬ 
larly large-bodied lycosid species or show be¬ 
havioural flexibility according to the available 
prey spectrum. Taking into account the need 
for sandy habitats for prey burrowing, A. rufa 
seems to be restricted to dynamic dune habi¬ 
tats in which soil development is poor and 
bare sand sufficiently available. Female A. fa¬ 
hrilis disappear from the population in late 
winter and early spring and are not available 
as host during the early beginning of the hunt¬ 
er wasp's adulthood (Bonte, unpub. data). 
Other large-bodied lycosid species, occurring 
in these habitats include only Alopecosa cu- 
neata (Clerck 1757) and A. harbipes (Sundev- 
all 1833) (Bonte et al. 2002), also having a 
spring activity, and are as a result neither 
available during summer. Possibly, A. rufa 
uses smaller lycosids as host during early 
summer (Arctosa perita (Latreiile 1799) or 
subadult A,. fahrilis ), potentially resulting in 
different spatial distribution patterns. 

The fact that A, rufa was only absent from 
small patches, inhabited by A. fahrilis seems 
to provide evidence that higher trophic levels 
require larger habitats (Holt 1996; Holt et al. 
1999), but may also be an artifact in our sam¬ 
pling strategy if smaller patches are only ac¬ 
cidentally visited by the spider wasp, hence 
reducing encounter chances in the field. Ac¬ 
cording to Holt (1999), higher trophic ranks 
need larger areas if dispersal is limited, have 
closed populations and show a high degree of 
specialization towards resource species. Cer¬ 
tainly the latter seems to hold for the (at least 
temporal and local) A. rufa-A . fahrilis asso¬ 


ciation. Assuming that our findings are not 
due to sampling artifacts, parasitoid absence 
in small populations indicates low dispersal 
abilities or at least low dispersal motivation in 
spider wasp within spatially structured suit¬ 
able habitat surrounded by hostile matrix veg¬ 
etation. 

Also, since both juvenile and adult A. fa¬ 
hrilis stay in their burrow during day, prey 
selection by A, rufa has probably to rely on 
olfactory and not on visual cues because bur¬ 
row entrances are not visible, as also observed 
in spider hunting wasps preying on day-active 
wolf spiders (e.g., Pompilus cinereus (Fabri- 
cius 1775); Bonte pers. obs.) and in sphecid 
wasps (although only at short distance after 
visual detection; Anton & Gnatzy 1998). 

In fragmented coastal dune habitats, host 
and parasitoid show, at least temporally and 
locally, a similar aggregative spatial popula¬ 
tion structure, indicating the importance of 
(lack of) dispersal and possibly spatially cor¬ 
related habitat characteristics in structuring 
patch occupancy patterns. Additionally, the 
spider hunting wasp A. rufa appears to behave 
as a specialized parasitoid during the main ac¬ 
tivity period of its optimal prey with restricted 
dispersal motivation. During this period, the 
presence of an assumed generalist parasitoid 
is a good predictor for the presence of noc¬ 
turnal and burrowing dune wolf spider in larg¬ 
er habitat patches in the Flemish coastal 
dunes. 
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